Theory of Fast Quantum Control of Exciton Dynamics 
in Semiconductor Quantum Dots 
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Optical techniques for the quantum control of the dynamics of multiexciton states in a semicon- 
ductor quantum dot are explored in theory. Composite bichromatic phase-locked pulses are shown 
to reduce the time of elementary quantum operations on excitons and biexcitons by an order of 
magnitude or more. Analytic and numerical methods of designing the pulse sequences are investi- 
gated. Fidelity of the operation is used to gauge its quality. A modified Quantum Fourier Transform 
algorithm is constructed with only Rabi rotations and is shown to reduce the number of operations. 
Application of the designed pulses to the algorithm is tested by a numerical simulation. 
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I. INTRODUCTION 

The possibility of controlling the dynamics of a quan- 
tum system has long captured the attention of workers 
in a wide range of physical systems. Quantum control 
can be realized by engineering a time-dependent Hamil- 
tonian which depends on a finite set of parameters. In 
quantum chemistry, this engineering has led to the pos- 
sibility of driving chemicaLreactions using tailored laser 
pulses and external fields.u Quantum control has been 
recently extended to semiconductor nanqatructures lead- 
ing, for instance- to controlled curraitsp coherent con- 
trol of excitonaJu and electron spin,El and controlled in- 
tersubband transitions of shallow donors using terahertz 
radiationJd 

In this paper we explore optical control of the ultimate 
quantum device in semiconductor nanotcchnology, i.e., a 
quantum dot. In a semiconductor quantum dot the elec- 
tronic levels have a density of states characteristic of a 
single atom. Yet, the dot is a mesoscopic system, i.e., in 
contrast to the single atom case, the quantization of the 
electronic levels is realized within a system that contains 
actually 10^ ^ 10^ atoms. A key ingredient in the quan- 
tum control of these semiconductor nanostructures is the 
robustness of the elementary excitation, the exciton. An 
electron-hole pair optically excited in an undoped quan- 
tum dot feels the presence of the large number of atoms 
in the material only through the static dieiectric constant 
and the electron and hole effective mass.Ll This allow us 
to treat excitons as excitations in giant atoms and to con- 
trol excitons with optical techniques similar to those used 
for the manipulation of atoms and molecules. However, 
unlike the atomic case, the dot is in a solid state envi- 
ronment, with the attendant decoherence. We shall also 
make use of the conduction band electron and the valence 
band hole as the constituents of an exciton. A quantum 
dot is like an empty box that can be filled with multiex-i 
citon complexes composed of many interacting excitons.u 
In these multiexciton states, the Coulomb correlation is 
taken into account and yet the spin configuration is trans- 
parent. The spin configuration can then be controlled by 



the light polarization of the optical pulses. 

The implementation of quantum algorithms is a par- 
ticular case of quantum control. The potentialities of 
semiconductor quantum dots in the implementation pJ 
quantum algorithms have been readily recognizecLt " Li 
as well as in conjunction with optical microcavities.E3~ti 
The use of optical control of excitons-,in dots for quan- 
tum operations has been suggested,E3 and a theory of 
the physical implementation of quantum algorithms ia. 
a dot using ultrafast optical pulses was investigated.E^I 
Ideas for a scalable quantum computer involving exci- 
tons in di|ff«:Giat dots and optical quantum control were 
proposed.li30'tl Advances in ultrafast optics in quan- 
tum dots make possible the manipulation of electronic 
excitations in a semiconductor nanostructure with time 
resolution in the femtosecond domain. So far, frequency 
selection is used to avoid unwanted transitions to states 
out of the computational space. Laser pulses of a narrow 
frequency range are too long in duration compared with 
the decoherence time for quantum operations. Thus, a 
design of fast control is necessary. A fast control allows 
us to make a reasonable number of operations well within 
the decoherence time. It may also be made an ingredient 
in the realization oLsophisticated error correction and de- 
coupling schemesJlj Wc will give an explicit design for the 
realization of fast control of two qubits encoded in two 
antiparallel-spin excitons in a single quantum dot. The 
slight increase in complexity of the optical setup is within 
the capability of the current experiments. An optimal 
design is an inverse problem to the finding a state given 
the Hamiltonian: the issue is to find a time-dependent 
optical electric field that produces a desired result in the 
shortest time as possible. The required experimental re- 
sources are realistic: lasers generating Gaussian pulses 
with two different frequencies that can be phase-locked. 
The synthesis of phase locked optical pulse fronv-sepa- 
rate femtosecond lasers has been recently reportedJ23 and 
here we propose an important application of this tech- 
nique. We explore the three different approaches to the 
control pipblem: an intuitive one making use of the area 
theoremEJ and also give an analytical tool, the cluster 
expansion of the evolution operation, much used in the 
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NMR spectroscopy,cHl and numerical optimization. 

The implementation of a 2-qubit Quantum Fourier 
Transform (QFT)E^ is used as a test case for the different 
methods of design. Any algorithm can be decomposed 



as a series of sin 
tional rotations '^^ 



e-qubit rotations and two-qubit condi- 
A series of two-color phase-locked opti- 
cal pulses is suggested to realize [the fast control of these 
fundamental rotations. FidelitjO is used to gauge the 
quality of the operations and of the complete algorithm 
of QFT within the decoherence time. Error correction 
could be added later to improve the result. 

The paper is organized as follow. Section |l| gives the 
structure of the multiexciton states from a microscopic 
model for a quantum dot. Section III explains the prin- 
ciples of three methods of design of optical pulses for a 
fast control in a subspace of multiexciton states and com- 
pares their results in fundamental quantum operations. 
Section IV contains a numerical simulation of the QFT 



algorithm in a quantum dot. The simulation takes into 
account the microscopic details of the laser-exciton dy- 
namics, including decoherence and the presence of mul- 
tiexciton levels outside the computational space. Details 
of the decomposition of the QFT in terms of only Rabi 
rotations for a general n qubit system are relegated to 
the appendix. Section ^ summarizes and draws a num- 
ber of conclusions. A brief description of the key idea of 
pulse shaping and the application to a different quantum 
algorithm were published in Ref. |l6|. 



II. MULTIEXCITON STATES 

The energies and wave functions of the multi-exciton 
states in a dot are calculated starting from two con- 
fined4evels of electrons and holes each in a parallelepiped 
QD.ca The electronic levels included are the first two 
states deriving from the localization of s-like conduc- 
tion band states. They carry a spin ±1/2. The hole 
levels derive from the localization of states in the p-like 
valence band heavy holes carrying a ±3/2 total spin in 
the direction of the growth axis. The size of the dot, 
40 nnix35 nrrui.5 nm, is typical of interface fluctuation 
quantum dotsEj Only Coulomb interaction between the 
carriers which conserves their conduction or valence band 
indices is taken into account exactly. This amounts to 
neglecting the electron-hole exchange which gives a fine 
structure of the excitonic levels depending on the sym- 
metry of the dot. We calculated this effect to be of the 
order of a few fJ,eV, which can thus be safely neglected 
in the discussion of fast control considered in this pa- 
per. Fig. 1^ shows the energy structure of the multiex- 
citon states. The multiexciton levels include zero, one, 
two, three, and four excitons in the dot. The choice of 
two levels each of electrons and holes limits the resultant 
number of excitons to four. The + or — refers to the po- 
larization of the light that has to be used to create each 
exciton. Only optically active multiexciton states are 



shown. Since the optically forbidden multiexciton states 
are not the source of unintended dynamics, they are re- 
moved from the following discussion. Our model adopts 
the measured dipole moment of 75 Debsie for the single 
exciton in a single GaAs fluctuation dotEZI and the tran- 
sition matrix elements between the multiexciton states 
are then calculated. The later values are used in the nu- 
merical simulations in section IV. Note that the values 



of the dipole moments in this kind of systems are one 
or two orders of magnitude higher than those of atoms. 
Theoretical estimates suggest that this giant dipole ef- 
fect seems to be stronger in quantum dots generated-by 
monolayer fluctuations than the self-assembled dots.E^ 
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FIG. 1. Energy levels of the multiexciton states in a quan- 
tum dot in meV. e+ — = 1764 meV. 



III. QUANTUM CONTROL OF MULTIEXCITON 
STATES 

A. Formulation of the Problem 

Quantum control consists in designing a time- 
dependent addition to the system Hamiltonian which will 
drive the state of the system from a prepared state to 
a designated state within a number of desirable condi- 
tions. In this paper, we shall focus on the dynamics of 
two excitons in a quantum dot. The controlling agent 
is a sequence of laser pulses. The dynamics of this two- 
exciton system not only serves as a powerful illustration 
of the more general case of multiple discrete states but 
also to form a set of "universal gates", i.e., fundamental 
operations in terms xif which any quantum computation 
may be constructedEj The system is not closed. In addi- 
tion to the laser interaction, the quantum dot has other 
excitonic states and its environment of the substrate and 
other dots is a source of decoherence. The time limi- 
tation due to decoherence and the resonance conditions 
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to avoid the unintended dynamics form two contradic- 
tory requirements under which the dynamics of the two 
excitons in a dot must be optimized. 

For the system of two excitons, we use a basis set of 
four states, in order, the ground state |0), the two ex- 
citons with the lowest energy at 1764 meV above the 
ground state and opposite polarizations denoted by |-|-) 
and |— ) and the biexciton state | — h) at 3527 meV. A 
a+ polarized photon can drive the excitonic transition 
|0) — > |-|-) and the biexcitonic transition |— ) \ — h). 
The two transitions have different energies due to the 
Coulomb binding energy of the biexciton. We can write 
the Hamiltonian of the four-level systems coupled to 
an external electromagnetic field with (t+ polarization, 
treated classically, in the form 



f^+(t)/2 
nX{t)/2 e+ 

e_ f^+it)/2 

fnut)/2 e_+ 



, (1) 



where = J2j ~ 'O) represents a time de- 

pendent Rabi energy provided by a train of phase-locked 
optical pulses. The dipole moment of the exciton |-|-) is 
denoted by d+ and / is a correction factor to the dipole 
moment in the exciton-biexciton transition matrix ele- 
ment due to Coulomb interaction. The amplitude of the 
electric field E+^j{t - Tj) = £+j{t - r,)e-*'^+(*-^^)e* 
is assumed to be slowly varying. As in the atomic case, 
the condition lu^ 3> rf+f+j- enables the rotating wave 
approximation used in 77+ above. Thus, the counter- 
rotating terms, such as Hq_ = ri!j_/2, are set to zero. 
Similarly the Hamiltonian associated with a (t_ polar- 
ized electric field is given by 



H- = 



n^{t)/2 

e+ fn^{t)/2 

n*_{t)/2 e_ 

fn*_{t)/2 £_+ 



(2) 



For simplicity of exposition we consider a sequence of 
non-overlapping pulses although in numerical simulations 
we have found possible to pack the pulses with 10% over- 
lap with negligible deterioration. Thus, we write the uni- 
tary time evolution operator from t = tq to t = T in the 
form 



UiT,To) 



N 



(3) 



where N indicates the number of pulses in the train and 
Tj,Tj^i the beginning and the end of the j-th pulse. For 
a given quantum operation [/(T, tq), the time optimiza- 
tion can be viewed as consisting of two components. The 
first is to have a minimum number of pulses N in Eq. (|^). 
Optical pulses can directly perform Rabi rotations with 
generators ax and ay but rotations with generator a^ 



need to be built as a combination of ay and ax ■ In our de- 
sign of the laser implementation of a quantum algorithm 
we try to decompose the required global transformation 
directly in rotations generated by ay and ax for both sin- 
gle qubit and conditional operations without appealing to 
Hadmard, C-NOT or conditional phase shift. We have 
demonstrated this by the construction of the Dcutsch- 
Josza algorithm££l and the Quantum Fourier Transform 
(see below). Since the saving is not exponential, in the- 
ory it may be considered trivial but in practice, especially 
in the initial stage of experimental implementation, the 
use of the right decomposition of the algorithms may be 
advantageous. 

The second component for a fast control is the time 
optimization of each pulse in the product of Eq. (|^), 
which is the main subject of this section. Consider the 
case of a fj-j- pulse. In the interaction representation, 
O = AO At denotes the transformed operator from O, 
with A{t) = e*'^"*, where Hq is a diagonal matrix with 
elements (0,e+,e_,e y). The term U^' in Eq. (||) be- 
comes for a+ pulse (with j understood below) 



(4) 



where V{t)'^+ is 



















(5) 



and A = e+ -|- e_ — e ^ is the biexciton binding energy. 

When only a circularly polarized light is used, Eq. (|^) 
shows that the four-level system behaves as a double two- 
level system, the first two-level transition (exciton tran- 
sition) being represented by |0) — > |+) and the second 
(biexciton transition) by |— ) — > | — h). 

Consider now the desired operation where the exciton 
transition is a Rabi rotation through angle a and the 
biexciton transition a Rabi rotation through a'. 



cos(a/2) 
sin(Q!/2) 







-sin(a/2) 
cos(a/2) 








cos(a72) 
sin(a72) 






-sin(a'/2) 
cos(a72) 



(6) 



The most direct solution for the realization of this trans- 
formation would be a two-pulse combination, 



E+{t) =£oe'^*/'^\-''^°+' 



(7) 



If the two pulses are resonant respectively with the two 

transitions, i.e. luq^ = e+ and coi^ = e ^ — e_, and 

sufficiently narrow in frequency, the pulse resonant with 
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the exciton transition would have neghgible effect on the 
biexciton transition and vice vpjjsa. However, this has 
been shown to be costly in tinie.li3 The problem is to find 
a composite pulse which would take much less time with 
tolerable deterioration of quality of the transformation. 

For the quality of the transformation, we follow Ref. |2j 
in defining the fidelity of the transformation as 



(8) 



where Ui is the ideal unitary operation, U is the uni- 
tary transformation generated by the optical pulses, and 
the overline denotes the average over all the possible ini- 
tial states. The operator WUi is denoted by / for short. 
The average over all the possible states is done by con- 
sidering an initial state with arbitrary complex coeffi- 
cients \ipin) = Cj\j) with the normalization constraint 
Si P — 1- The fidelity can be then written in the form 



ijkl 



(9) 



and, in the four-level system considered here, the over- 
line average is then on a hypersphere S in deter- 
mined by the normalization condition. This average 
{l/S) Jg (Pci(Pc2(Pc3(fc4C*Cjclci is easily evaluated in 
polar coordinates and gives 



F = 1/10 + 1/20 ^(/../;^. 



(10) 



The difference of the coefficients from those of Ref. ^ is 
due to their additional restrictions on the coefficients Cj . 
Our choice gives a more conservative estimation of the 
error in the operations. 



2. The average Hamiltonian method 

The cumulant expansion (also known as the Magnus 
expansionnj) of the evolution operator C/J^ in Eq. (^) is 
given bycj 



(12) 



The first term of the expansion corresponds to a time 
average of the interaction Hamiltonian, 



Vi = I dtV{t). 
'o 



The second term is given by 



V2^— dt dt'[V{t),V{t')] 



(13) 



(14) 



Keeping only the first term in the exponent constitutes 
the average Hamiltonian approximation. An estimation 
of the error in the truncation of the cumulant expansion 
is given by the second term. 



3. Numerical approach 

The parameters in Eq. (^ are varied to find the max- 
imum fidelity. To lessen the numerical effort, physical 
considerations guide the reduction of the number of pa- 
rameters varied. The first two approximation methods 
are also useful as starting points. 



C. Examples of pulse design 



B. Pulse Design 

In this subsection, we explain three different ap- 
proaches to pulse design to shorten the time of the quan- 
tum operation. 



1. Approximation by the Area Theorem 

In the limit of very long pulses, the area theoremill 
determines the intensity of a Gaussian pulse that has to 
be used for a given rotation a 



s^/nd^ 



(11) 



For a single two-level system the pulse width s in Eq. (|1] 
can be made arbitrarily small, but in the four-level case 
we are strongly limited by the resonance condition to 
1/s, £od+ <C A. In order to shorten the time duration of 
the whole pulse, an intuitive approach would be to allow 
the two components of Eq. (^) to overlap in frequency 
but keep each satisfying the area theorem. 



We illustrate the above methods for a single qubit 
operation, i.e., a parallel rotation of both the exciton 
and biexciton transitions. For simplicity, let / = 1 and 
si = s. Both theoretical estimates and experimental 
measurements have the / value not far from unity. In 
any case, the extension to / 7^ 1 can be made in a simi- 
lar manner to the treatment on the conditional rotation 
given below. We consider a composite pulse by superpos- 
ing and phase-locking the the two pulses in Eq. (^ with 
£0 = £1 a-nd ujQ^ = e+ and uji^ = e+ — A. It remains 
to choose a value for £q{s) by each of the three methods 
above and tests its efficacy by evaluating the fidelity of 
the operation. 

In Fig. 11(a) the fidelity for a = a' = tt rotation is 
plotted as a function of the temporal width of the Gaus- 
sian pulse s. The corresponding value for the peak of the 
Rabi energy fig = d+£o{s) is also given in Fig. ||(b). The 
value of the biexcitonic binding energy A is 1 meV. The 
results by the area theorem approximation are shown as 
the dashed lines. The fidelity is close to unity only for 
s ^ 1/A, corresponding to a region where the frequency 
selectivity is preserved. If for instance a 98% Fidelity 
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is required, the area theorem approach will lead to op- 
tical pulses with s > 4 ps. The area theorem is not the 
best procedure of time optimization for single-qubit op- 
erations. 




Area Theorem 
Numerical Maximization 
Average iHamiltonian 



0.2 



s (ps) 



FIG. 2. (a) Fidelity as a function of the temporal width 
of the Gaussians s for a parallel rotation of a = a' = tt. 
(b) Peak value of the Rabi energy Qq = d+£o{s)- Dashed 
lines: the area theorem approximation. Dotted lines: the 
average Hamiltonian approximation. Solid Lines: numerical 
maximization of the fidelity. 

Applying the average Hamiltonian approximation to 
the restricted pulse specified above leads to the single- 
qubit rotation C/'^+ in the form of Eq. with chosen 
values for a = a' and for s, leading to So given by 



sV^d+(l + e-(^«/2)^) 



(15) 



The Gaussian term in the denominator on the right gives 
a correction to the area theorem, Eq. ([Tl]). The results 
are shown as dotted lines in Fig. ^(b). An estimate of the 
error of the average Hamiltonian approximation may be 
made by evaluating the second-order term in the cluster 
expansion given by Eq. (14). A rough estimate is pro- 
vided by replacing the Gaussians with square pulses of 
width s, 



V2 ^ (d+£o/A)2(sin(As) - As 

2Ascos(As/2) - 4sin(As/2))H 



(16) 



where H is a diagonal matrix with ele- 
ments Expanding in the limit of short 

pulses As < 1 we get (d+£o/A)2(As) V3. The 

correction to the area theorem in the first order term, 
Eq. (|l^), is by contrast ~ (d+£o/A)(As). Faster pulses 
make the lowest order Vi ^ ¥2- The resultant fidelity 
by the average Hamiltonian method is shown as dotted 
lines in Fig. ||(a). Note that it is possible to obtain a 98% 
Fidelity using much shorter pulses, of the order of 100 fs. 
In the limit of very short pulses this correspond to pulses 
spectrally very broad which do not distinguish between 
the two transitions but yield a nearly parallel rotation. 



The results of the numerical maximization using one 
variable Eq by Brent's methodtJ are plotted as solid lines. 
The optimal curve £q{s) deviates considerably at short 
times from the area theorem approximation but is close 
to the average Hamiltonian approximation throughout 
the whole range of s. 

The second example is a conditional operation for two 
qubits, viz., a (t+ biexcitonic transition without affecting 



|0), i.e. 

— TT. 



a rotation Uj in Eq. 



the excitonic \+) 

with a ~ Q and a' — tt. For the combined pulse in 
Eq. (|7|) we consider now (p — n, and again £0 = Si and 
^0+ = = e+ - A. 

From the average Hamiltonian approximation (the first 
order term in the cluster expansion), we obtain relations 
for the three parameters of the pulse £q, s and si for the 
desired rotations. 



a = rf+£oV^(s-sie"('^"i/2) )^ 
a' = d+£o%/^(si-se-('^^/2)'). 



(17) 
(18) 



For a given value of si, the other two parameters may be 
solved in the case with a = and a' = tt, 



s = Sle-(^^^/2)^ 
£0 = V^/rf+(si-se-(^^/2)^ 



(19) 
(20) 



In the limit of large A the solution gives s ^ elim- 
inating the term resonant with the excitonic transition 
and £0 ^ \^/sid^ in accord with the area theorem 
for the biexcitonic transition. For A ^ this sys- 
tem has always a solution for any a a' . Correction 
to the average Hamiltonian approximation may be es- 
timated in analogy to the parallel rotation case in the 
limit si, s ^ 1/A and give for V2 a diagonal matrix with 
elements (-^, ^, where 

-^1 - + 2(A^i)' - 3(As)(AsO^], 



32' A 

h - ^i^nAs.f + 2iAsf - 3(Asi)(As)2]. 



In Fig. H we show (a) the fidelity and (b) the peak 
Rabi energy for the a = Q and a' ~ tt transformation, 
for all three methods. The area theorem approximation 
amounts to taking a single pulse resonant with the biexci- 
ton transition. For the numerical maximization we max- 
imize the fidelity for a given si value as a fuaction of 
s and £0 using the downhill simplex method.E2l We see 
clearly that the average Hamiltonian again gives a very 
good approximation: the deviations from the numerical 
maximization are negligible in most of the region. Also in 
this case we see that the use of a composite pulse provide 
a considerable saving in the time for the operation. 
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FIG. 3. (a) Fidelity as a function of the temporal width si 
of the biexciton Gaussian component in the composite pulse 
for a rotation of a' = tt only for the biexciton transition, 
(b) Peak value of the Rabi energy. Dashed lines: the area 
theorem approximation with a single pulse resonant with the 
biexciton transition. Dotted lines: the averaged Hamiltonian 
approximation. Solid Lines: numerical maximization of the 
fidelity. 
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FIG. 4. (a) Fidelity as a function of the temporal width s\ 
of the square pulse and of the biexciton Gaussian component 
of the shaped pulse for a rotation of a' = vr only for the biex- 
citon transition, (b) Peak value of the Rabi energy. Dashed 
lines: when a single square pulse resonant with the biexci- 
ton transition and the area theorem is used. Dotted lines: 
composite Gaussian pulse using the average Hamiltonian ap- 
proach, same as the dotted lines in Fig. H. 



As a last example, we investigate vifhether a single 
square pulse shape can serve the function of the two over- 
lapping pulses. For square pulse an exact analytical ep. 
pression for the C/"*" can be given. It has been suggestedEil 
that off-resonant unwanted transitions can be corrected 
using square pulses. In the specific case discussed above 
this corresponds to the use of a single pulse resonant with 
the biexciton transition with s\ satisfying the conditions 

siA = V4m2 - X-kU , (21) 
sil7+ = TT , (22) 

with integer m. Eq. (|2^ ) gives a tt rotation for the biexci- 
tonic transition in accord with the area theorem, and the 
condition in Eq. ( |2l| ) sets to zero the off-diagonal terms 
in the 2x2 block corresponding to the excitonic transi- 
tion. However, additional phases in the diagonal cor- 
responding to a rotations for the exciton transitions 
are introduced which decrease the fidelity of the opera- 
tion. We calculate the fidelity and peak Rabi energy for 
a conditional tt rotation using a single square pulse res- 
onant with the biexcitonic transition as functions of the 
temporal width si of the square pulse and compare it 
with the shaped pulse result of the average Hamiltonian 
approximation in Fig. ^. In the comparison, note that 
si in the square wave is the temporal width but in the 
shaped pulse is the half width of the biexciton Gaussian 
component. The fidelity of the square wave shows oscil- 
lations with maxima roughly corresponding to the con- 
ditions in Eqs. ( |2l| ) and (^) but never reaches as high 
as the two-pulse case. Moreover, the spread in frequency 
of the square pulse spectrum is a source for unintended 
dynamics for higher exciton energy levels in the physical 
dot, while Gaussian pulses avoid this problem. 



IV. QUANTUM FOURIER TRANSFORM 

The theory of control of the two excitons will now be 
applied to construct a physical implementation of a quan- 
tum algorithm, the two-qubit Quantum Fourier Trans- 
form. One qubit is given by the presence or absence of a 
polarized exciton in the dot and the second qubit by 
the a~ exciton. For the single-qubit operation we have 
to act both in the exciton and in the biexciton transi- 
tion. The case of parallel tt rotation discussed in the 
previous section corresponds, therefore, to a single-qubit 
operation on the first qubit. As a conditional two-qubit 
gate we use a C-ROT operation, which is essentially the 
C-NOT operation with a tt rotation replacing the logical 
NOT operation. The C-ROT is the conditional dynam- 
ics of adding to the quantum dot a (t+ exciton only if an 
exciton with (t_ polarization is ahead v, there. This corre- 
sponds to a (7+ biexcitonic transitional without affecting 
the excitonic |-|-) — > |0). The second example discussed 
in the previous section is a C-ROT controlled by second 
qubit i.e. by the ct_ exciton. Single qubit and conditional 
rotations are easily generalized to arbitrary angles. The 
exact mapping for the four-level system to two qubits is 
given by 

{|0),|+>,|-),|-+)}^{|00),|01),|10),|11)} (23) 

Unlike the NMR implementations, it is not possible here 
to make use of the free evolution of the interacting qubits 
since it requires keeping track of the oscillation at the 
optical frequency and, therefore, an optical control over 
several picoseconds with sub-femtosecond accuracy. By 
working in the interaction representation, we get rid of 
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this drift term, making the design of the control more 
convenient. At the end of the sequence of pulses for a 
given algorithm the interaction representation does not 
affect the computation since the read-out is always done 
in a eigenstate of the system. Therefore the control of 
the qubit is always active and is constructed in terms 
of rotations with ax and ay generators between pairs of 
levels in the four-level system. 

The Quantum Fourier Transform (QFT) is the key in- 
gredient to a number of importaiit.quantum algorithms, 
in particular Shor's factorizationHj Consider an n-qubit 
state 



■ xi), where x — Xi2^ 



(24) 



The QFT is defined as a linear operator on an orthonor- 
mal basis of |0), . . . , - 1), where N = 2", with the 
action 



* q=Q 



(25) 



analogous to the construction of Bloch states in a one 
dimensional chain. The standard way to construct QFT 
employs two basic operations: the Hadamard gate on 
the jth qubit Hj and the conditional phase gate 
where j is the control qubit and fc is the target qubit .E3 
The two-qubit QFT can be realized for instance by the 
simple sequence H^S^^iHi- However, if we decompose 
each of the three operations in Rabi rotations we end up 
in using more optical pulses than necessary. In fact, each 
Hadamard transformation requires at least two optical 
pulses, i?j(7r, i:)i?j(7r/2, y), where i?j(0, e) the rotation on 
j-th qubit in the e direction with angle Q. Following the 
decomposition in Ref. 32, wc find 5i2 with the sequence. 



i?i(-|,y)Ci?or2,i(f,x)Ci?OT2,i(-f,i)i?i(|,f) 

i?i(f,y)i?2(-f ,y)i?2(f ,i)i?2(f ,y) . (26) 

In CROTj k, j is the control qubit and k is the target. 
For example. 



CR0T2,i{e,x) = 



10 

10 

cos{e/2) -isin(6l/2) 

.0 -isin(6'/2) cos(6'/2) 



(27) 



and the bar over suffix 2 indicates a rotation on the target 
only for the control qubit in the state 0, 



CROT2 ii9,x) 



cos(6i/2) -isin(6i/2) 

~ism{e/2) cos(6'/2) 

10 

1 



The total number of pulses for the QFT is then 12. 
We redefine the QFT as 

Umqft — BUqft'S, (29) 
where S is the all-qubit inversion (xi — > 1 — Xi), 

n 

E\x) = \x), where X ^"^{1- Xi)2'~-\ (30) 



and B is the transformatipa (xi 
be termed boustrophedon,c^ 



B\Xr. 



■X2X1) 



\X1X2 ■ 



), which may 



(31) 



(28) 



In the Appendix, we prove that Umqft is a composi- 
tion of rotations of generators a^ and ay for states of 
any number of qubits, denoted by Umqft- By avoid- 
ing the pulse-consuming Sij, this saves time by using a 
smaller number of pulses than Uqft- Umqft can be 
used directly in phase estimation or factorization algo- 
rithms without the need for B and E, the global qubit 
transformations which are just re-labcling of the qubits. 
In a physical implementations there is the possibility to 
make global qubit transformations that are simple rela- 
beling, at no cost from the point of view of the quantum 
control. If for instance a quantum computer is composed 
of a chain of one-half spins, at any time we can decide to 
fiip all spin up in spin down and vice versa. This all-bit 
inversion is a simple relabeling. We do not need to apply 
any pulse to the chain; we have just to remember that 
in the readout. The same can be done by switching in 
reading the string of qubits from the right to the left in- 
stead than from the left to the right, which corresponds 
to the boustrophedon transformation in Eq. (^l|). Al- 
though this saving in time is of the order polynomial in 
n, for the current attempt at physical implementation of 
prototype quantum computers it could provide a helpful 
simplification of the experimental procedure. 

For n = 2, TV = 4 the pulse sequence for Umqft is 

CROT2A\.m2{-\.v)R2{^.mi{-\.v) ■ (32) 

We carried out a numerical simulation of the dynamics 
of the multiexciton levels for this MQFT algorithm with 
and without the use of composite pulses. We took the 
peak of the Rabi energy to be 2 meV, larger than the 
1 meV binding energy of the biexciton. The width of the 
pulses is calculated using the area theorem approxima- 
tion and the average Hamiltonian. The corresponding 
values of fidehty for Umqft are 0.257 and 0.992. The 
pulse sequence is completed within 6 ps. 

In order to check the robustness of the use of composite 
pulses in the presence of dephasing, we include the spon- 
taneous emission in the simulation by adding the Lind- 
blad operators in the equation of motion for the density 
matrixEj 
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where, 



4 

E 



L,pL]--L]L,p--L 



(33) 



Li = Vf|0)(+|,L2 = Vf|0)(-|, 
i3 = Vf|+)(- + |,L4 = \/f|-)(- + |, (34) 

r = being chosen to approximate the measured 

dephasing time.E3 These operators represent all the pos- 
sible spontaneous emission pathway in the four-level sys- 
tem. 

There are many equivalent ways to solve the master 
equation in terms of a nonlinear stochastic differential 
equation for a normalized state vector We choose to 
use the quantum state diffusion (QSD) equation^ 



\di}) ^ --H\i,)dt 



^(i,-(L,))|^)d77, 



{L]){L,)]\r^)dt 



(35) 



where {L) = {tl)\L\tl)) and rjj are independent complex 
random variables. The density matrix can be expressed 
as p = M\'ip){ip\ where M denotes ensemble average 
and the expectation value of any operator O is given by 
M{ip\0\'4>). Inclusion of dephasing in this way reduces 
the fidelity for the shaped pulse sequence of MQFT from 
0.992 to 0.892. 



V. CONCLUSIONS 

In the quantum control of multiexciton states in semi- 
conductor quantum dots, we have shown that the use of 
composite pulses makes possible the realization of quan- 
tum operations in time scales of the order of a hundred 
femtoseconds. In addition to the theory of methods of 
constructing the pulses, we gave explicit examples to help 
experimental implementation. We adopted the concept 
of fidelity as a measure of the quality of a pulse sequence. 
We showed how to construct a sequence of pulse based 
only on the physical (7^ and ay rotations. A numerical 
simulation of the application of the shaped pulses to the 
two-qubit Quantum Fourier Transform in a single semi- 
conductor quantum dot provided a test of the pulse shap- 
ing. While the work so far provides a complete blueprint 
for an experimental demonstration of a simple quantum 
computation, future work for a more realistic computer 
includes the inter-dot for scaling up the system, design 
of optical control to minimize decoherence, and design of 
optical implementation of quantum error corrections for 
digital control of decoherence and unintended dynamics. 
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APPENDIX A: 



In this appendix the pulse sequences for the MQFT are constructed for an arbitrary number of qubits. From 
Eqs. (p5|) and d29), the action of MQFT and the inverse are given by, 



N-l 

UmqftIx) = ^ E e'"'^'/^l9>, 



(Al) 



9=0 

N-l 

uIqftIi) = ^ = E e-^^^'^^/^k). (A2) 



Define Ry to be a y-rotation on all the qubits. Then, 



AT-l 



Ry\x)^l[R,i-n/2,y)\x) = ^ E e^.- ^^''^■^^"^^■^b). (A3) 



Now consider the combined transformation 

N-l N-l N-l 

RyUlcFTk) - E Ry\^){^\Ul,QFT\<i) = ^ E E e^U-^'^'^-^^^e-'-'^^/^lp) (A4) 

x=0 p=0 a:=0 



P— J — 1 a^j— 

= E n e-"*«"'/^ (cos(7rg(^)/2)<5,^,p^ + z sin(7rg(^'V2)<5,„i-p,) b) (A6) 

p=0 j=l 

N-l n n 

= E n e"''' n (cos(7rg(-')/2)J,^,p^ + zsin(7rg(^")/2)'5,„i_p,) b) (A7) 
= E e'"^' n {^'"''^'^^^^ + e""'=^''5<?.,o) n (cos(7rg(^)/2)<5,^,,^ + z sin(7rg(-'-)/2)<5,^.i-p,) b>- (A8) 

p=0 fc=2 j=l 

where we use the definition g^^' = J2k=j+i The ak is defined through the relation YVj=i e~'^***^' = nfe=2 

and a = X]fc=2 "^fe- Note that afe depends only on k and N. 

The two products in the last line ( [A^ ) may be related to the rotations, 

n M 

ni?.K,^)k> = n {^'""^'3,^,1 + e-^'"'^'K,o) \q), (A9) 

i=2 j=2 

n n 

Rj{-Trq^''>,x)^Rj{-TT E gfc2J■-^i)= II CROTk^ji-nV-^x). (AlO) 
These relations lead via 



n n 



^f^l/QFT = e'"/2 n n CROTk,j{-nV-\x)\[Rk{ak,z), (All) 

j=l k=j + l fe=2 



to the conclusion that 



n n 



C/AfQFT = e-^"/' n n CROTk,,{nV-^x)Y[Re{-ae,z)Y[ Rmi-7r/2,y). (A12) 

j=l /c=j+l e=2 m=l 

Note that CROTkj may be moved to the right as in the circuit diagram Fig. |^ but not past any rotation involving 
the target qubit j. Finally, by using R{a, z) = i?(— 7r/2, y)R{a, x)R{tt/2, y) we obtain a pulse sequence which involves 
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only rotations and conditional rotations in the x, y direction to implement MQFT. Note that the total number of 
operations for n qubits is O(n^). 



R. 



R. 



R, 



R, 



R. 



R, 



R, 



R. 



RzH 



B 



FIG. 5. Circuit diagram for QFT, Eq. (A12), with the operations in the order from left to right. Each horizontal line 
represents a qubit. The operations are explained in the text. The ones connecting two quibit lines represent logic gates of 
controlled rotations. 
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